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Abstract 

We use higher-dimensional bosonization to study the normal state of elec- 
trons in weakly coupled metallic chains interacting with long-range Coulomb 
forces. Particular attention is paid to the crossover between Luttinger and 
Fermi liquid behavior as the interchain hopping t± is varied. Although in 
the physically interesting case of finite but small t± the quasi-particle residue 
does not vanish, the single-particle Green's function exhibits the signature of 
Luttinger liquid behavior (i.e. anomalous scaling and spin-charge separation) 
in a large intermediate parameter regime. Using realistic parameters, we find 
that the scaling behavior in this regime is characterized by an anomalous 
dimension of the order of unity, as suggested by recent experiments on quasi- 
one-dimensional conductors. Our calculation gives also new insights into the 
approximations inherent in higher-dimensional bosonization; in particular, we 
show that the replacement of a curved Fermi surface by a finite number M 
of flat patches can give rise to unphysical nesting singularities in the single- 
particle Green's function, which disappear only in the limit M — * oo or if 
curvature effects are included. We also compare our approach with other 
methods. This work extents our recent letter [Phys. Rev. Lett. 74, 2997 
(1995).]. 
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I. INTRODUCTION 



The properties of the normal metallic state of correlated electrons in d = 1 dimension 
are usually summarized under the name Luttinger liquid behavior The single-particle 
Green's function of a Luttinger liquid exhibits several striking differences to the Green's 
function of a conventional Fermi liquid: the absence of a coherent quasi-particle peak, spin- 
charge separation, and anomalous scaling properties characterized by interaction-dependent 
power laws. Non-Fermi liquid behavior has recently been observed in several experiments. 
First of all, the normal-state properties of the high-temperature superconductors cannot be 
interpreted in terms of a conventional Fermi liquid picture [|| . Motivated by this experimen- 
tal fact, Anderson H proposed that even in d = 2 the single-particle Green's functions of 
correlated electrons can exhibit Luttinger liquid behavior. This unconventional normal-state 
behavior was then shown to be a possible key to understand the novel superconducting state. 
In fact, the interlayer-tunneling theory of high-temperature superconductivity advanced by 
Chakravarty, Anderson and co-workers |4j] is based on the fundamental assumption that in 
the normal state the single-particle Green's function G(k, to) satisfies for wave- vectors k suf- 
ficiently close to the Fermi surface and for sufficiently low frequencies to (measured relative 
to the chemical potential) an anomalous scaling law of the form 

G{k a + sq,su) = s^ 1 G{k a + ci,uj) , (1) 

where k a is a wave-vector on the Fermi surface, and 7 > is some interaction-dependent 
exponent, the anomalous dimension. A finite value of 7 is one of the fundamental charac- 
teristics of Luttinger liquid behavior, whereas in a Fermi liquid 7 = 0. It should be stressed, 
however, that a generally accepted microscopic derivation of Eq. (0) for interacting fermions 
in d > 1 does not exist. 

Another class of materials where non-Fermi liquid behavior appears to have been ob- 
served experimentally consists of weakly coupled metallic chains [[J, which are based on 
highly anisotropic conductors. Although at sufficiently low temperatures T these systems 
have the tendency to develop various types of long-range order, above the ordering tem- 
perature the metallic state displays clearly non-Fermi liquid behavior. The interpretation 
of recent photoemission studies of these systems in terms of the Luttinger liquid picture 
leads to values for the anomalous dimension 7 in the range 1.0 ± 0.2 13-0] . Because in 
one-dimensional lattice models with short-range interactions (such as the Hubbard model) 
the anomalous dimension is always small compared with unity ||], the large experimentally 
measured value of 7 suggests that the long-range nature of the Coulomb interaction must 
play an important role in these systems. Of course, the experimental systems are not strictly 
one-dimensional, because the interchain hopping t± is small but finite, and even for t± = 
the electrons on different chains interact with three-dimensional Coulomb forces. We thus 
arrive at the problem of weakly coupled metallic chains, which is the central topic of this 
paper. 

In the last years the problem of coupled Luttinger liquids has been studied by many au- 
thors [|T0| |23|] . Given the fact that in the absence of interchain hopping the non-perturbative 



bosonization approach can be used to calculate the single-particle Green's function in a 
controlled way, some authors [|12],^3|,|17|,^0|,^2],^3| have attempted to supplement the one- 



dimensional non-perturbative bosonization solution by some kind of perturbation theory 



2 



in powers of t±. One disadvantage of this strategy is that in this way interchain hopping 
and electron-electron interactions between the chains are not taken into account on equal 
footing. In this work we shall therefore adopt a different strategy: Instead of combining 
one-dimensional bosonization with perturbation theory in powers of t±, we shall use the 
recently developed higher- dimensional generalization of the bosonization method [f24]-|3"0f to 
treat interchain hopping and interactions between the chains on equal footing. Thus, our 
approach is non-perturbative in t± and in the interaction. In particular, for t± = we 
recover the well-known one-dimensional bosonization result for the Green's function of the 
Tomonaga-Luttinger model. 

The rest of this paper is organized as follows: In Sec. [TT] we briefly describe the higher- 
dimensional bosonization approach. In particular, we emphasize the approximations inher- 
ent in this approach and discuss its limitations. In Sec. [TTI| we shall show how the problem 
of electrons that are confined to one-dimensional chains (i.e. without interchain hopping) 
and interact with three-dimensional Coulomb forces can be solved in a straightforward way 
with higher-dimensional bosonization. Of course, this problem can also be solved by means 
of the conventional one-dimensional bosonization technique |3T]-|3~3|[ . However, as shown in 
Sec. [IV], our approach can also handle the case of finite interchain hopping t± without having 
to rely on an expansion in powers of t± or in the interaction. Although the quasi-particle 
residue is found to be finite for any t± ^ 0, we show by explicit calculation that there exists 
a large intermediate regime where the equal time Green's function satisfies an anomalous 
scaling relation similar to a Luttinger liquid. In Sec. [V] we compare our results with those 
obtained by other approaches. Finally, in Sec. |V| we present our conclusions. 



II. HIGHER-DIMENSIONAL BOSONIZATION 

In this section we shall give a brief summary of the higher- dimensional bosonization 
approach. The foundations of this novel non-perturbative approach to the fermionic many- 



body problem go back to ideas due to Luther p3 and Haldane |25| . More recently, a number 



of authors have applied this technique to problems of physical interest [[26H 29I1- The basis of 
higher-dimensional bosonization is the subdivision of the Fermi surface into a finite number 
of patches Pg, a — 1, ... , M. To be specific, let us consider in Fig. [I] a typical Fermi surface 
associated with a periodic array of chains coupled by weak interchain hopping. A possible 
subdivision into M = 8 patches is shown in Fig. Each patch P£ is then extended into 
a 3-dimensional squat box [^BJ as shown in Fig. |3[ The characteristic size A of the patches 
should be chosen sufficiently small so that within a given patch the curvature of the Fermi 
surface can be locally ignored. On the other hand, the cutoffs A and A cannot be taken 
arbitrarily small, because higher- dimensional bosonization becomes only useful in practice if 
it is possible to ignore momentum transfer between different patches (the so-called around- 
the-corner processes [23]). Of course, this is only possible if the interaction is dominated 
by forward scattering. More quantitatively, we assume that the two-body interaction / q 
becomes negligibly small if |q| is larger than some cutoff q c <C min{A, A}. The geometry is 
shown in Fig. || Assuming that the nature of the interaction is such that the above condition 
can be satisfied, it is reasonable to make the following two approximations [B8LB9|: First of 



all, we ignore momentum transfer between different boxes (diagonal-patch approximation). 
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The relative number of matrix elements neglected in this case is in d dimensions of order 
gf/(A d_1 A) <C 1. The second fundamental approximation is the local linearization of the 
energy dispersion. Measuring wave- vectors locally with respect to coordinate systems located 
at points k Q in the centers of the patches (see Fig. §), a general energy dispersion ek may 
be expanded as 

(Qu) 2 (a a ) 2 

where q° = qv a and q" = q— (q-v a )v a . Here v Q = v Q /|v a | is a unit vector in the direction 
of v a . Note that in general the masses mjj* and m" depend on the patch- index. Because in the 
grand-canonical formulation of statistical mechanics the energy dispersion appears only in 
the combination e^ — fi, the chemical potential \i cancels the constant in Eq. (g) provided 
k a is chosen to lie exactly on the Fermi surface. We now ignore all terms in Eq. (|j) that are 
quadratic and higher order in q, i.e. we approximate £° ~ v a • q. In particular, the term 
(q°) 2 /(2m") is neglected. Note that retaining this term would lead to patches with non- 
zero curvature. Once we accept the validity of the above approximations, the single-particle 
Green's function can be calculated without further approximation in arbitrary dimension. 



As shown in Refs. Vm and 291, the result for the Matsubara Green's function at finite 



inverse temperature (3 can be written as 



G(k a + q,iid n ) = J dr J dre-^-^G^r,^^'^ , (3) 



where 



1 „i(q-r-o) n r) 

^^ w ^-—^ , (4) 

and the Debye- Waller factor Q a (r, r) is given by 

Q a (r,r)=R a -S a (r,T) , (5) 

with 

oat \ 1 V- f RPA (^ iuJ rn) cos(q • r - u m r) 

{,) ~WL W m - V • q) 2 • (?) 

Here tb n = 2n(n + \)/ (3 are fermionic Matsubara frequencies, and u) m = 2nm//3 are bosonic 
ones. f RPA (ci, iu) m ) is the screened interaction in random-phase approximation (RPA), which 
is given in terms of the bare interaction / q and the non- interacting polarization IIo(q, iw m ) 
via the usual relation 

/^(q,^ m )= /q , . (8) 

1 + / q n (q, iu) m ) 
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In the RPA screened interaction Eq. (§) only the long wavelength limit of Ilo(q, icJ m ) enters, 
consistent with the neglect of the around-the-corner processes. Thus, the q sums in Eqs. 
(f§)-(0) are implicitly restricted to the regime \q\\ \ £ A, \q±\ £ A. However, as long as the 
external wave-vector q in Eq. (|j) is sufficiently small, we may ignore these cutoffs. Then it is 
easy to see that Gq(t, t) in Eq. (^) is proportional to §( d ~ 1 \r±), where the d — 1-dimensional 
vector r_i_ consists of the components of r that are perpendicular to v a . In fact, for V and 
(3 — ► oo the integration in Eq. (PS) is easily performed analytically, with the result 

where rn = r • v". Because of the prefactor 5 < - d_1 - ) (rj_), we may replace r — > riiV Q in Eqs. ([?[) 
and (|^). Furthermore the r-integral in Eq. (|J) is effectively a one-dimensional one, because 
the integrations over the components of r perpendicular to v a can be done trivially due to 
the 5-function in Eq. (^). Of course, this simplification is an artefact of the linearization. 
The modifications of the above results due to the non-linear terms in the energy dispersion 
have recently been calculated in Ref. |3(J. Most importantly, for non-linear energy dispersion 
one should replace in Eqs. (f|) and (0) 

I I 



[lOJr, 



(10) 



where £° is given in Eq. (0). Note that the non-linear terms remove the double pole in Eqs. 
(P) and (|7|). Moreover, non-linearities in the energy dispersion lead also to the replacement 
of the function Gq(t,t) in Eq. (|J) by an interaction-dependent Green's function Gf(r, r), 
which does not involve the singular prefactor 5^ d ~ 1 \rj_) given in Eq. For an explicit 
expression of G"(r,r) see Refs. and 



The above non-perturbative expression for the single-particle Green's function are valid 
in arbitrary dimensions and for arbitrarily shaped Fermi surfaces. Furthermore, in d = 1 
these expressions correctly reduce to the well-known bosonization result of the Tomonaga- 
Luttinger model. We now apply the above results to study the problem of coupled Luttinger 
liquids. 



III. THE COULOMB-INTERACTION IN CHAINS WITHOUT INTERCHAIN 

HOPPING 

Before discussing finite interchain hopping t±, it is instructive to study first the case 
t± = 0, where the electrons are confined to the individual chains. Of course, electrons 
on different chains still interact with the three-dimensional Coulomb interaction, so that 



this is not a purely one-dimensional problem. The latter has been discussed in Ref. |p2 
Thus, making the continuum approximation for motion parallel to the chains, the Fourier 
transform of the bare interaction is given by 

U = e 2 al dr x J2^— , (11) 
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where the r^-sum is over the two dimensional lattice of chains, with interchain lattice spacing 
a±. The prime indicates that the = term must be properly regularized (see below). 
For |q| <C aj 1 Eq. (|TT| ) reduces to the familiar result / q = 47re 2 /q 2 . Note that we do not 
replace by hand the long-range Coulomb interaction by an effective screened short-range 
interaction. The screening problem will be solved explicitly via our bosonization approach. 

Because for t± = the particle number is conserved on each chain, the problem can be 
solved by means of standard one- dimensional bosonization techniques pT|-p5|. However, the 



solution can also be obtained quite elegantly within the framework of higher-dimensional 
bosonization. In the absence of interchain hopping, the Fermi surface consists of two parallel 
planes, as shown in Fig. |. These planes can be identified with the patches discussed in Sec. 
U Let us label the right plane by a = +, and the left one by a = — . The (linearized) energy 
dispersion is simply given by £ q = avFq x , where v F = |v a | is the Fermi velocity. As usual, a 
dimensionless measure for the strength of the interaction is F q = z// q , where v = 2/(7rt>^a 2 L ) 
is the density of states at the Fermi energy (the factor of 2 is due to the spin-degeneracy). 
An important length scale in the problem is set by the Thomas-Fermi screening wave- vector 
k, which can be defined by rewriting the dimensionless interaction F q in the long wavelength 
limit as _F q = K 2 /q 2 . This yields k = a^y/Wng, where the dimensionless coupling constant 
g is given by 



e 2 



9 = • (12) 

hvf 

To evaluate Eqs. ©-(0), we need to know the RPA interaction in Eq. (§), which involves 
the non-interacting polarization n (q, iu> m ). In the absence of interchain hopping the po- 
larization has the usual one-dimensional form, which is in the long wavelength limit given 
by 



n (q, iw m ) = v - — ~ — = 7^ — rr? « • ( 13 ) 

2 • q- iuj m 1 + UJ^I \v F q x y 



For better comparison with the calculations for finite interchain hopping presented in Sec. 
|TV| , it is convenient to express the RPA interaction in terms of the dynamic structure factor 
in the usual way (we take the limit (3 — ► oo) 

f RPA {^ iu m ) = / q - fl r duS RPA {^ u)-^ , (14) 

"'0 u z + uj^ 

with 

1 ^ vr \l + / q n (q, W + *0+)J 1 ; 

Using Eq. (|13D we obtain for to > 

S rpa {^u) = Z (1 5{lo-uj ci ) , (16) 
with the collective mode and the residue given by 
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1 + FqV F \q x \ 
w F \q x \ 



2 A /l + F q 7raiJl + F q 



(17) 
(18) 



Substituting these results into Eq. (H) it is then easy to show that 



which for V 



R° 

oo reduces to 
1 



--Yf 



R a = -- 



(w q + |v« • ql) 1 



F 2 
q 



o q x \2 A /l + F ( 



(19) 



(20) 



BZ 



Similarly, we obtain 



ReS a (x,r) = -- 



dq x 
2 Jo q x 



cos(q x x) 



1 + Ea 



l + F a 



y/l+i r q«Ffe|T| 



ImS Q (x,r) 



sgn(rj 







sin^x) 



«F9se T 



— e 



-v F q x \r\ 



BZ 



(21) 



(22) 



where x = ar x , and where for any function h(q) the symbol (h(q)) BZ denotes averaging 
over the first transverse Brillouin zone (BZ), 



(Kq)) 



1 



BZ 



I 2tt \2 



dq y / a± - dq z h(q) 



(23) 



To show that the g^-integrals in Eqs. ( ^Tp and (^) exist at large q x without the addition 
of an ultraviolet cutoff it is necessary to know the behavior of _F q for large |q|. Trivially in 
the approximation F q = K 2 /q 2 the interaction falls off as |q| -2 . If one takes into account 
that the one-dimensional continuous and two-dimensional discrete Fourier transform ( [TT]) 
of the potential shows the same behavior (see below) it is easy to show from Eqs. (|2~T|) and 
(p2|) that after the BZ averaging the integrands for r = fall off like (n/\q x \) 5 . Because the 
(^-integrals exist, k plays the role of an ultraviolet cutoff. For finite r the integrands in Eqs. 
and ( |22|) only fall off like (n/q x ) 2 but the g x -integrals still exists. 
Let us first consider the case r = 0. Then we need to calculate the following Brillouin 



zone average 



F 2 
q 



2\2Jl + F q L/l + F q + l] 



BZ 



1 + £k 
l + Fa 



BZ 



(24) 



In the regime g <C 1 the Thomas-Fermi screening length k 1 is large compared with the 
transverse lattice spacing a±. Because in this case all wave- vector integrals are dominated 
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by the regime |q| £ k it is allowed to use the continuum approximation / q = 4-7re 2 /q 2 for 
the Fourier transform of the Coulomb potential. The averaging over the transverse Brillouin 
zone in Eq. can then be done analytically, with the result 



Hence, 



icb(qx 



Q a (x,0) 



2itvf 



7<d|/« + ^1 + (q x /K) 2 



(25) 



2-KVp Jo 



dq x 



1 - cos(q x x) 



(26) 



For x 3> n 1 it is now easy to show that 



Q a (x, 0) ~ -j cb \n[Kx] 



(27) 



with the anomalous dimension given by 



7 cb = lim n jcb(q x ) 
qx-*0 



2-KVi 



(28) 



The logarithmic growth of the static Debye- Waller factor is one of the characteristics of a 
Luttinger liquid. Note that in a strictly one-dimensional model the long-range Coulomb 

m 



interaction leads to a Wigner crystal phase 



and not to a Luttinger liquid. In this 



case the anomalous dimension diverges in the thermodynamic limit. However, the Coulomb 
interaction between electrons on different chains removes this divergence, and leads to a 
Luttinger liquid |32| with modified spectral properties for k ^ kp (see below). 

We would like to point out that Eq. fl2"B| ) is only valid for g «C 1, where 7^ -C 1. It would 
be incorrect to extrapolate this result to the regime where g is of the order of unity, which 
is experimentally relevant. In order to calculate the anomalous dimension in this regime, 
the Fourier transform / q of the interaction has to be properly calculated. For larger g we 
have to take in to account the lattice structure in the transverse direction. To calculate 
the one-dimensional continuous and two-dimensional discrete Fourier transform ( |TT| ) of the 
Coulomb potential it is necessary to regularize the = 0-contribution as in the strictly 
one-dimensional case ||34|| . The characteristic feature of the regularized = 0-term is its 
logarithmic divergence — cln(a g x ) for small q x , where ao is the typical extension of the 
one-particle wave-function in the transverse direction and c is a constant. The prefactor 
and the behavior at larger momenta depend on the special regularization procedure chosen. 
This divergence is responsible for the divergence of the anomalous dimension in the strictly 
one-dimensional model with Coulomb interaction, which leads to the Wigner crystal phase 
]34| of the one-dimensional electron gas. For g C 1 we already saw that the chains coupled 
by the three-dimensional Coulomb interaction are not in the Wigner crystal but in the 
Luttinger liquid phase. In the following we show that for arbitrary g the term — cln (aoq x ) 
for an infinite number of chains is exactly cancelled by a term cln(a±q x ) from the rj_ 7^ 
part of the Fourier transform (0). Therefore the system is a Luttinger liquid for all g. 
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In our regularization procedure we assume that the one-particle wave function in the 
transverse direction is given by a Gaussian distribution corresponding to a harmonic con- 
finement potential. We replace the r_|_ = 0-component in Eq. (|TT| ) by 

Q 2 6^ f°° f°° g — *9a: T "a;g~ r x/ a 



dr_L / 

J — 



a 7r J-oo J-co [ r £ -f r^j ) 



a 2 e 2 
ag7r 



f°° dr ± K (\r ± \q x )e~ r y a o, (29) 

J—oo 



where Kq(x) denotes the modified Bessel function. The logarithmic divergence of Kq for 
x — > leads to the logarithmic behavior for small g x discussed above. 

To obtain a fast convergence of the Fourier transform flTT| ) for rj_ ^ 0, we use the Ewald 
summation technique |[35|| . With the help of the Theta function transformation 



V^ e -Qi|r x l 2 -iqxTx _ ^ ^ c -|qx-Gx| 2 /(4a)^ 
r ± a^CY 

where Gj_ denotes a vector of the reciprocal lattice in the transverse direction and the 
integral 

4= rd ge -w^iv= , 1 , (3i) 

0T./o J r 2 + | r± |2 



one obtains 



exp{- [q 2 x + 


q±-G ± 


2 ]/(4«o)} 




q±-G ± 


2 



G 

da 



+ 2a\e 2 r — E e- a l r ^l 2 -^^ 

Jan a _ j n 



2a 2 ± a 2 e 2 P° daa \— e -^ 2 ) 

'o 1 + ao a 



+2aie 2 r^_i_ e -d/(^). (32) 
Jao a 1 + agcr 

For a G [0, oiq] we used in Eq. ([H|) the transformation Eq. fl3~0D from the direct to the 
reciprocal lattice. Now cto can be chosen in such a way that in the sums in the first and 
second line of Eq. (^) only a few lattice vectors contribute. In the following we always 
choose ao = 4/ao. For this cto the sum in the second line in Eq. (J53) can be completely 
neglected and in the first line only the vectors Gi = (m y 2ir/a±,m z 2TT/a±) with m y and 
m z € [—5, 5] contribute. The last two terms in Eq. fl3"2| ) contain the rj_ = 0-term of Eq. 
(PI). From Eq. ([!2|) it is easy to see that only for |q| <C a J 1 is is allowed to approximate 
the Fourier transform of the Coulomb potential by its continuum limit 47re 2 /q 2 . 

To determine the anomalous dimension from Eqs. ( p4|) and (|28|), we have to take the 
q x — > limit in Eq. (|32f) . For q x = we obtain 
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0,q ± )=47re 2 ]T 

G, 



+2a 2 , e 2 



2a 2 t e 2 In 



exp { — 


|q± - 


Gi 


lV(4«o)} 




|q± - 


Gi 


12 


f°° da 

Ja a 


£ e 


— a|r 


L | 2 -iq ± -R ± 



4^ 



(33) 



i.e. no logarithmic divergence. As already discussed above, the logarithmic divergence of 
the r_|_ = 0-term is exactly cancelled. Therefore the anomalous dimension is finite and can 
be calculated numerically following Eqs. (f2_|) and (|28|). In the regularization procedure we 
introduced the new parameter clq. Physically it is clear that only ratios a±/ao > 1 are 
relevant for quasi-one- dimensional systems. In Fig. [7] we present numerical results for ^ c b 
with a±/aQ = 3 (dotted line) and a±/ao = 6 (dashed line). For comparison we also show 
7 c fc using the continuum approximation / q = 47re 2 /q 2 (solid line). From Eq. ( |32|) it is easy 
to show, that F q falls off as K 2 /q 2 for large |q|. Therefore our above discussion concerning 
the existence of the integrals Eqs. fl2~l|) and (|22|) is also applicable to the lattice Fourier 
transform. 

In a Luttinger liquid the momentum integrated spectral function p(uj), which apart 
from a one-electron dipole matrix element determines angular integrated photoemission, is 
algebraically suppressed near the chemical potential |36| 



p(u) oc M 7c \ (34) 

Recent photoemission studies of quasi-one-dimensional conductors suggest values for the 
anomalous dimension in the range 1.0 ± 0.2 This would be hard to reconcile with 

a model involving short range interactions, as e.g. the anomalous dimension in the one- 
dimensional Hubbard model never exceeds 1/8. A rough estimate of the dimensionless 
coupling g for experimental relevant parameters following Ref. || leads to values of g of the 
order of one. As seen from Fig. [7] our treatment using realistic Coulomb forces quite easily 
leads to anomalous dimensions in the experimental range. 

In the following we discuss the momentum dependent spectral function p < (q x ,u). In the 
strictly one- dimensional model with an interaction potential which is finite at q = it is 
known that for anomalous dimensions 7 smaller than 1/2 a power law divergence occurs at 
the energy vpq and for 7 < 1 at v c q, where v c is the long wavelength limit of the g-dependent 
charge excitation velocity (spin-charge separation) [^7J,^]. From Eq. (JT7J) one obtains for 
the coupled chains a q^-dependence of the charge velocity 



u c (qj.) = Km -5 = Vp Ji + F q± . (35) 
qx ^ u Qx 

It is not obvious from Eqs. ([H]) and (p2|) whether a charge peak occurs or not, because the 
charge-velocity t> c (q_i_) is subject to the BZ averaging and the singularity can be washed 
out. On the other hand, because the spin velocity Vp is independent of q^, one obtains for 
7ch < I a sharp threshold singularity at uj = —vp\q x \, just like in one dimension, 

p%(q*,u) ~ e{v F [q x - k F ] - uj){v F [q x - k F ] - uj)^~ 1 ' 2 . (36) 
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For 7 c b > | the singularity is washed out, but the threshold survives. On the Fermi surface 
the BZ averaging is irrelevant and one obtains for 7^ < 1 



p<(A; F ,cu)~e(-^)(-^- 1 (37) 

as in the one- dimensional Tomonaga-Luttinger model. In order to evaluate the momentum 
dependent spectral function for arbitrary q x and lu, we have to use numerical techniques 
developed in Ref. [|39|] . It turns out that a broadened charge peak occurs. For increasing 
anomalous dimensions, i.e. increasing g, the charge peak broadens. To clearly demonstrate 
spin-charge separation, we present in Fig. 8 spectra for 7 cb = 0.33, corresponding to g = 0.8, 
with a_|_/a = 3. The spectra are calculated for large but finite systems. The system size 
is given by the dimensionless level spacing v K = (2it/L)/ k. In Fig. |8| we chose v K = 1/80. 
For technical reasons we multiplied the Green's function by a factor exp {— x\t\} before we 
performed the time-frequency Fourier transformation. This leads to a Lorentzian broadening 
of the discrete spectral function and smears out the threshold behavior at u) — Vp{k x kp). 
As can be seen from Fig. [5], the artificial broadening x/{ v fk) = 0.01 is much smaller than 
the natural broadening of the charge peak. The power law behavior of the spectra at 
uj = vp(k x — kp) leads to an asymmetry of the related peaks in Fig. [| We checked that for 
a given \ t ne numerical spectra do not change if the system size is further increased. In the 
numerical calculations we have to cutoff the g x -integrals in Eqs. (|2l| ) and (^) at q x = nn, 
with some number n. In Fig. 8 we chose n = 5 but we checked that for the momenta and 
energies presented the spectra are the same for larger values of n. We found that the charge 
peak for —0.5 < {k x — kp) / k < disperses like uj = v c {k x — kp) with an effective g dependent 
charge velocity v c which has to be determined numerically. For anomalous dimensions larger 
than one no divergence at u = vp{k x — kp) occurs and the spectra are dominated by a very 
broad charge peak. 

To summarize this section, we showed that the system of one-dimensional chains coupled 
by the long ranged Coulomb interaction is a Luttinger liquid with an anomalous dimensions 
of the order of one. Furthermore the charge singularity of the strictly one-dimensional 
Luttinger liquid is broadened due to the BZ averaging, but is still present. In the following 
we extend the model and include transverse interchain hopping. 



IV. FINITE INTERCHAIN HOPPING 

A typical Fermi surface of an array of chains with small finite interchain hopping is 
shown in Fig. [TJ. In the quasi-one- dimensional materials discussed in Ref. [0 the intrachain 
hopping t x is an order of magnitude larger than the interchain hopping t y = t±, which is 
again a factor of 10 larger than t z . We therefore assume transverse hopping only in the y 
direction. In this section we shall apply our higher- dimensional bosonization approach to 
this problem. We would like to emphasize that we are ultimately interested in realistic Fermi 
surfaces without nesting symmetries. Unfortunately, by replacing a curved Fermi surface of 
the type shown in Fig. [I] by a finite number of locally flat patches (see Fig. 0), we introduce 
unphysical nesting symmetries. In order to obtain physical results that can be compared 
to experiments, the singularities caused by these artificial symmetries should therefore be 
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separated from the physical plasmon mode in the dynamic structure factor. In this section 
we shall show how this can be done in practice. 



A. The M-patch model 

In the higher- dimensional bosonization approach with linearized energy dispersion the 
Fermi surface is approximated by a finite number M of flat patches P£. In order not to 
break the inversion symmetry of the Fermi surface, we choose an even number of patches 
of equal size such that for each P£ there exists another patch P£ such that the local Fermi 
velocities satisfy v a = — v a (see Fig. for M — 8). Then the non-interacting polarization 
n (q, z) is at long wavelengths given by 

2u M/2 (w a ■ a) 2 

= i^o^ • (38) 

where it is understood that the sums are over all patches with v* • q > 0. For finite M the 
poles of the RPA interaction, which are simply the zeros of 

l + / q n (q,z) = , (39) 



can be easily obtained by plotting the right-hand side of Eq. ( |3"8D as function of uj 2 where 
uj = z is real, and looking for the intersections with — l// q . For generic q all (v° • q) 2 are 
different and positive, and we can order < (v ai • q) 2 < (v a2 • q) 2 < . . . < (v aju / 2 ■ q) 2 . A 
repulsive interaction / q leads then to zeros (u; 2 )( a ), a — 1, . . . , M/2 of Eq. (|39| ) as function 
of z 2 lying between the unperturbed poles, 

(v ai • q) 2 < (cu 2 )« < (v" 2 • q) 2 < . . . < (v Q «/ 2 • q) 2 < {u\f M ^ . (40) 

In the large M-limit the first M/2 — 1 solutions form the particle- hole quasi-continuum, 
while the largest mode w*^ 2 can be identified with the collective plasmon mode. For q- 
values which lead to v ai • q = and (or) [y ai ■ q) 2 = (v" J • q) 2 the number of poles in the 
quasi-continuum with non-vanishing residue is reduced. Denoting by Z q the residue of the 
plasmon mode, the dynamic structure factor is given by 

S RPA (q,u J ) = Z^(u J -uJ+Sgr o t(q,u) , (41) 

where S^^(q, uj) is due to the contribution of the other modes, which merge for M — > oo 
into the particle-hole continuum. In the absence of interchain hopping ^^(q, a;) = 0. 
Furthermore, the plasmon mode and the associated residue are given in Eqs. (O) and (15). 



1. The plasmon mode 

The crucial observation is now that for small but finite t± and M —>■ oo the low-energy 
behavior of the Green's function is still dominated by the plasmon mode, which therefore 
describes the crossover between Luttinger and Fermi liquid behavior. Thus, to leading order 
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in 9, we may ignore the contribution S^^(q,u>) in Eq. (|4lD. In this case the constant part 
R a of the Debye- Waller factor is still given by Eq. (19), but now with tu q and Z q given by 
the plasmon mode of the M-patch model. In the strong coupling limit it is easy to obtain 
an analytic expression for u; q and Z q . Anticipating that at strong coupling there exists a 
pole ojq of Eq. 
term is 



]) with uj 1 = 0(F q ), we may expand Ilo(q, z) in powers of z . The leading 



n (q, z) 



2u 



M/2 



M ^ 



a=l 



v ■ q 



1 + 



V • q 



+ 



q 



Substituting this approximation into Eq. 



it is easy to show that 



2 M/2 



\2 _|_ Eq= ! ' 

a=l \ 



M/2 V*. q ) 4 



E 



q) 



+ O (1/F q ) 



and 



2VF, 



9 A//2 

— V(v c 

a=l 



qj 



1/2 



Introducing these expressions in Eq. (p!9|) we obtain to leading order in the coupling 



ql + rj 2 q 2 ± 5 [(jJ — Vf 



(42) 



(43) 



(44) 



(45) 



where 



M/2 



M 



(46) 



Of=l 



is a small parameter for a weakly corrugated Fermi surface corresponding to small interchain 
hopping. Note that q^ in Eq.(f45D denotes the q y - and g 2 -components of q, and should not 
be confused with the vector q" in Eq.(Q). The contribution to the constant part of the 
Debye- Waller factor is given by 



pa 




(47) 



BZ 



To leading order the patch dependence drops out. For r] — 0, corresponding to zero interchain 
hopping the integral is logarithmically divergent. For finite rj the q x integration obtains a 
finite lower cutoff which leads to a contribution In [«;/ (77 1 qx | ) ] - To leading order in In (1/ 77) 
one can take the limit q x — > in F q . This yields 

1 



K pl ~ 



In (I/77) lim 
2 q x ^0 



BZ 



(48) 



The prefactor of the logarithm is just the strong coupling expression of the anomalous 
dimension of the zero hopping 77 = limit. Now, if is finite, the Riemann-Lebesgue 
Lemma [|(| guarantees that Spi(r\\\ a , 0) goes to zero for increasing ry. Therefore the system 
is no longer a Luttinger liquid, but a Fermi liquid with a finite quasi-particle weight P8]f2T[] 

Z a = e Ra . (49) 
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2. The nesting mode 



Clearly, by approximating the realistic curved Fermi surface shown in Fig. [I] by a col- 
lection of flat patches, we have introduced an artificial nesting symmetry. For example, the 
patches PI and P^ +4 (with a = 1,2,3,4) in Fig. | can be connected by constant vectors 
which can be attached to an arbitrary point on the patches. Simple perturbation theory 



41 1 indicates that this nesting symmetry gives rise to logarithmic singularities, leading to 



a breakdown of the Fermi liquid state. However, unless there exists a real physical nesting 
symmetry in the problem, this singularity has been artificially generated by approximating 
a curved Fermi surface by flat patches. 

Let us first show how these nesting singularities manifest themselves within our bosoniza- 
tion approach, and then give a simple quantitative argument how these singularities disap- 
pear due to curvature effects. The crucial observation is that for a given q there exists one 
special patch Pf such that the energy |v^ • q| is much smaller than all the other energies 
|v a ■ q|, a (3. Anticipating that for sufficiently small q\\ = ir 13 ■ q there exists a 5-function 
peak in the dynamic structure factor at the nesting mode cx |v^ ■ q|, we see that for 
small q\\ the energy is much smaller than all energies |v a • q| with a ^ (3. Hence, the 
energy dispersion of the nesting mode can be approximately calculated by setting z 2 = in 
all terms with a ^ (3 in the sum Eq. (|3lf). This yields in the regime of wave- vectors defined 
by 



\v p ■ q| < |v a • q| for all a ^ (3 



(50) 



the following expression for the polarization 



n (q, z) 



V 

M 



M-2 + 



2(v^ ■ q 



V ■ q) 2 



(51) 



The collective mode Eq. ( |3~9"D is then easily solved, with the result that the dispersion of the 
nesting mode is given by 



l + F a 



\1 + 



M-2 rp 
M r 1 



l v/3 ■ q| 



(52) 



For the associated residue we obtain 



v\vP ■ q| 



M 



1 J_ M ~ 2 TP 

1 + Hw-^q 



(53) 



Using Eq. 



we see that the nesting mode gives rise to the following contribution to the 



constant part of the Debye- Waller factor associated with patch Pjf , 



R 



nest 



— £' 



F 2 
q 



1 + 



M-2 
M J 



[1 + ^ 



1+F a 



i M-2 £T 



+ 1 



2 ' 



(54) 
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where the prime on the sum indicates that the sum is over the wave-vector regime defined 



in Eq. (50). To understand the geometric meaning of Eq. (50), let us note that for our 
M-patch model the angle between v' 3 and the neighboring velocity v^ +1 is of the order of 
0/M, where 
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M 



(55) 



Setting q\\ = • q and q± = ■ q (where is a unit vector orthogonal to ir 13 ) it is easy 
to see that the condition Eq. (|50|) is equivalent with -C 6\q±\/M. We conclude that for 
V — * oo the integral in Eq. fl54D is proportional to 



R 13 oc — 

±l/ nest ^ jyj- 



dq_i_ 



Q i±l M dq\\ 



(56) 



which is infrared divergent. Obviously the logarithmic divergence is removed if we take the 
limit M — ► oo. It should be kept in mind, however, that in the limit M — > oo the size of the 
patches vanishes, so that the condition q c <C min{A, A} is violated. As discussed in Sec. [T| 
this condition is necessary to justify the neglect of the around-the-corner processes in the 
higher-dimensional bosonization approach, see Fig. |3|. Nevertheless, because for vanishing 
patch cutoff we recover a curved Fermi surface, the above calculation suggests that the 
nesting singularity will disappear as soon as the finite curvature of the patches is taken 
into account. In Sec. |1V B 2j we shall use the higher-dimensional bosonization result for 
curved patches [p9j , j3Cfl to show that this is indeed the case. Here we would like to give 
a simple quantitative argument which captures the essential physics. Suppose we retain 
the quadratic terms in the expansion Eq. (^]) of the energy dispersion close to patch P^. 
Obviously the term qj_/(2m±) (describing the curvature of patch Pf) becomes important 
for Uf|<7||| ~ q'j_/(2m±), where for simplicity we have assumed that m x = m±_ and |v^| = Vp 
are independent of the patch index (3. Hence, we expect that for curved patches the lower 

2 

limit for the (^-integral in Eq. flSB]) will be effectively replaced by 2 \m^\ VF ■ We conclude that 
the effect of curvature can be qualitatively taken into account by substituting 



r0g±/M dq\\ 
Jo q\\ 



e i±/ M dq\\ 
«j — = m 

2\m±\v F "II 



2\ui±\vf0 
q±M 



(57) 



In physically relevant cases we expect \m±\ ~ kp/ivpO), so that the right-hand side of Eq. 
(|57|) reduces to the integrable factor of ln[2A;_p/(M|gj_|)]. A more rigorous justification for 
the regularization given in Eq. (|57|) will be given in Sec. [IV B 2. 



B. The 4-patch model 

We now confirm the general results discussed in the previous subsection by explicit 
calculations for M = 4, where the collective mode equation is quadratic and can be solved 
exactly. Let us start with the Fermi surface shown in Fig. [|. The local Fermi velocities are 
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(e^ cos 6 + e y sin 9)vf 
(e x cos — e y sin 9)vf 
(—e x cos 9 + e y sin 
(— e x cos — e y sin 6')vj? 



and the non-interacting polarization is now 



n i v " ^ (v° ■ q) 2 

" ,W ^£ (v. q )H< 



The collective mode Eq. (|39|) reduces to the bi-quadratic equation 



i + ^)(^ + S)^ 2 + (i + F q )^ = o 



where we have introduced the notation 

£ q = v 1 • q = Vp(q x cos 9 + q y sm( 
£ q = v 2 ■ q = Vp(q x cos 9 — q y sin( 



Eq. (BUI) is easily solved, 



2 



+ P2 
Sq T Sq 

2 

q \ Sq 1 Sq 



+ 



+ (l + F q )(^-4 



2\2 



2\2 



1/2 



1/2 



The right-hand side of Eq. (|64|) is non-negative because 



1 | gAg+j. 
2 / 2 



1 

4 



q | 2 



F 2 ( ^_^, +( i+F q )(^-a 2 



(l+^q«q>0 . 



(58) 



(59) 



(60) 



(61) 
(62) 



(63) 



(64) 



(65) 



Both modes co> q and co> q give rise to ^-function peaks in the dynamic structure factor. We 
obtain for uo > 



S RP (q, Cfc>) = Zq5(^ - CJq) + Z n 5(uj - CJq) 

with the residues given by 



,/ , ,2^±|| C2t2 

y ^q~ ^q 

2 ^q ^q - ^q 

^ SqSq ^q 2 

2ci} r , U! 2 — u 2 

q 



(66) 

(67) 
(68) 
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In the limit 9 — > we have £ q — > £ q — > wfQ'x, so that co> q — > yl + -Fqi^l^l and u; q — > u.f|9x|- 
It is also easy to see that the residue Z q in this limit reduces to the result given in Eq. (jl8l) , 
while the residue Z q vanishes. To examine whether interchain hopping destroys the Luttinger 
liquid behavior, it is sufficient to calculate the static Debye- Waller factor. Combining Eqs. 
(H), (0) and flSfiP, it is easy to show that 



g a (r,|V a , 0) = [l - cos(v Q ■ qr,,; 

q 



cos(v Q • qry) /J 



(cj+|v a -q|) 2 



+ 



v a -q|) 2 (£ q +|v a -q|) 2 



(69) 



1. The plasmon mode 

Let us first focus on the first term in the second line of Eq. (|69|). Because for t± — > 
this term smoothly reduces to the corresponding expression in the absence of interchain 
hopping, the mode u; q can be identified with the physical plasmon mode discussed in Sec. 
[IV A| . From Sec. |TTT] we know that the Debye- Waller factor is essentially determined by the 
long wavelength regime |q| & k. In this regime F q ^> 1, so that we may use the strong- 
coupling approximation for u; q and Z q . From Eqs. (|BTf ) and it is easy to see that, up to 
higher orders in 9/F^, the collective density mode u; q can in this regime be approximated 
by 



^ q ^ FV /l + F qV /g 2 + ^ 2 , . (70) 
Substituting this expression into Eq. (|67|), we obtain 



■W«s + *«i (71) 



2,/l + F„ 



Comparing Eq. (|7T| ) with Eq. flT%l), we see that the only effect of the interchain hopping 
is the replacement \q x \ — > ^Jq% + 9 2 q 2 . Thus, the plasmon mode c<j q yields the following 
contribution to the constant part of Eq. (|69"D for small 9 




*3 = -o /. <M ; q • (72) 



If we set 9 = we recover the previous result in the absence of interchain hopping (Eq. 
(P0|)), which is logarithmically divergent. This divergence is due to the fact that for 9 = 
the first factor in Eq. (|72|) can be pulled out of the averaging bracket. However, for any 
finite 9 the q x - and (^-integrations are coupled, so that it is not possible to factorize the 
integrations. Hence, any non-zero value of 9 couples the phase space of the q-integration. 
Because for 9 — > the integral in Eq. (|T2|) is logarithmically divergent, the coefficient of 
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the leading logarithmic term can be extracted by ignoring the (^-dependence of the second 
factor in the averaging symbol of Eq. (|72|) . Then we obtain to leading logarithmic order 



R a , ~ — 
pl 2 Jo 



dq x 



"2 k 



lim 

q x — >0 



= hm 

y 



F 2 
q 



2J1 



BZ 



F 2 
q 



2J1 + F q /1 + F q + 1 



BZ 



-Jcb 



hi 



+ &1 



(73) 



where 7 c fe is given in Eq. (|28[), and &i is a numerical constant of the order of unity. 

The contribution of the plasmon mode to the spatially varying part S a (r»v a , 0) of the 
Debye- Waller factor at equal times can be calculated analogously. Note that r\\ = v a ■ r = 
±r x cos 9 ± r y sin 9. Repeating the steps leading to Eq. ([72]), we obtain 



^y(niv a ,o) 



1 

2 Jo 



x cos(g x r|| 



cos 



ql + 2 ,/i + F a ,/i + F a + 1 



BZ 



(74) 



Because the Thomas-Fermi wave-vector k acts as an effective ultraviolet cutoff, the value 
of the integral in Eq. ( |74D is determined by the regime |q| ^ k. For < 1 we may 

approximate in this regime cos^^ry) w 1 under the integral sign. Furthermore, for ft|r||| ^> 
1 the oscillating factor cos(g x r||) effectively replaces k by |r||| _1 as the relevant ultraviolet 
cutoff. We conclude that in the parametrically large intermediate regime 

(75) 



(76) 



K 1 < r||| 






logarithmic order 






^(niv a ,o)~- 7cfe 




+ b 2 



where b 2 is another numerical constant. 



2. The nesting mode 



From the general analysis presented in Sec. [IV A| we expect that the second mode o) q in 
Eq. ( |69| ) will give rise to an unphysical logarithmic growth of Q a (r\\v a , 0) at large distances, 
which is caused by the nesting symmetry of the Fermi surface in our simple 4-patch model. 
We now verify this expectation and then use the results of Refs. [^9|,[3(J to show by explicit 
calculation that curvature effects remove the nesting singularities. For convenience we choose 
the integration variables q\\ = v 1 • q = q x cos 9 + q y sin 9 and q± = —q x sin 9 + q y cos 9. Then 
£q = vpq\\ and £ q = Vp(q\\ —29q±) to leading order in 9. Note that the condition |g||| ^ 9\q±\ 
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is equivalent with |£ q | £ |£ q |. Geometrically this means that the wave- vector q is almost 
parallel to the surface of the first and fourth patch, so that its projection v 1 • q = — v 4 ■ q 
on the local normals is much smaller than the projection on the normals v 2 and v 3 of the 
other two patches. In this regime we obtain from Eq. ( |64T) to leading order 



and from Eq. (I 



l + F. 



f7 v f\Q\\\ 



q \1 + 

^ ^j-|g||| 
^4[1 + ^]I[1 + F q ]^ 



kill ~ o\q±\ 



We conclude that for \q»\ & 0\q±_\ 



(77) 



(78) 



v 



(w q +|v«-q|) s 



4^| g|| |[l + f]f[l + F q ]l[ I^l 



1+^ 



l] s 



(79) 



Note that the above expressions agree with Eqs. Q52|)— fl53|) if we set there M = 4. The 
factor of l/\q\\\ in Eq. ([f9| ) implies that the static Debye- Waller factor grows logarithmically 
for rii — > oo. However, for realistic Fermi surfaces of the form shown in Fig. |I| the nesting 
symmetry responsible for this behavior does not exist. To remove this artificial nesting 
symmetry, let us now replace the completely flat patches of Fig. [5] by the slightly curved 
patches shown in Fig. || The corresponding energy dispersions can be taken to be £° = 

„a 

v" • q + with negative effective mass m±. For weakly coupled chains we should choose 
|m_i_| pa m\\/0, where m\\ = kp/vF is the effective mass for motion along the chains. As 
shown in Refs. |29| , [30fl , one important effect of the curvature terms on the Debye- Waller 
factor is the replacement given in Eq. flTUD , which removes the double pole in Eqs. ([|) and 
([?]). Taking this effect into account, we see that the contribution of the nesting mode to the 
constant part R a of Eq. (B9l) becomes 



r: 



v q 



2sgn(e) 



[oo. 



+ I£q1) 



(80) 



with Z q and o) q given in Eqs. (68) and (0). From the above discussion we know that 
possible nesting singularities are due to the regime |g||| ^ 0\q±\. Thus, restricting the limits 
of the integral in Eq. ( |30|) to this regime, we obtain from Eqs. ( j64|) and (|68|) in the strong 
coupling limit 



f 2 a z q 



VF\q\\\ 
y/2v F \q\\\ 



(81) 
(82) 



Recall that for the three-dimensional Coulomb interaction the strong-coupling condition 
ufq ^> 1 is equivalent with |q| <C k. Putting everything together, we find that the contri- 
bution from the critical regime |g||| & 0\q±\ to Eq. (|50"D can be written as 
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_2 

y/2K f K , ImJ r e \i±\ , |?[||sgn(?|| — 



„2 

The (^-integration can now be performed analytically. The integral is proportional to — — 



which cancels the sing ular factor of ^ in Eq. (g3|). We obtain 



K "' ~ ~ y/2(y/2 + 1)WVV F L ^ 



In ( 2| "^ M ) + 6 3 



|mx| ' 



(84) 



where 63 is a numerical constant of the order of unity. This is the same type of integral as 
in Eq. fl57|), so that our simple intuitive arguments given at the end of Sec. |IV A| are now 



put on a more solid basis. As already mentioned, in physically relevant cases we expect 
\mx_\vp6 « kp, so that we finally obtain Rn est ~ — lebbi where 64 is another numerical 
constant of the order of unity. Thus, for patches with finite curvature the contribution of 
the nesting mode is finite. It is also easy to see that the curvature terms do not modify the 
logarithmic small-6 1 behavior of R a due to the plasmon mode given in Eq. fl73[). This is so 
because the leading ln(l/#)-term in Eq. ( [73]) is generated by the energy scale Vp9\q±\, which 
is by assumption larger than the curvature energy q]_/{2\m±\). Comparing our result for 
Rnest with the contribution from the plasmon mode given in Eq. ([75]) , we see that for small 
9 the contribution of the nesting mode does not modify the leading logarithmic behavior 
Rp t ~ — 7 c fe ln(l/#) due to the plasmon mode. Hence, to leading logarithmic order, the static 
Debye- Waller factor is dominated by the contribution from the plasmon mode, so that we 
may write R a ~ RZ, and similarly for S a (r\\v a , 0). 

3. Anomalous scaling 

Because R a is finite for any non-zero 6, the system is a Fermi liquid, with quasi-particle 
residue 

Z a = e Ra oc 9^ b . (85) 

Thus, for 9 — > the quasi-particle residue vanishes with a non-universal power of 9, which 
can be identified with the anomalous dimension of the corresponding Luttinger liquid that 
would exist for 9 = at the same value of the dimensionless coupling constant g. Combining 
Eqs. (|73| ) and (|76|), we obtain for the total static Debye- Waller factor 



Q a (r\\v a , 0) = R a - S a {r\\v a , 0) 



m(«|r|||) + &l , k- 1 < |r,|| < (9k)~ 1 , (86) 



where b is a numerical constant of the order of unity. Exponentiating this expression, we 
see that the interacting Green's function satisfies the anomalous scaling relation, 

G a (r/s,0) = s 3+7c "G a (r,0) , k~ x < |rn | , |r,| \/s < (9k)- 1 . (87) 
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In momentum space this implies for u = and 6 k <C |q| C k an anomalous scaling law 
of the form given in Eq. ([!]) with 7 = 7^. Thus, in spite of the fact that the system is 
a Fermi liquid, there exists for small 6 a parametrically large intermediate regime where 
the interacting Green's function satisfies the anomalous scaling law of Luttinger liquids. 
Moreover, the effective anomalous exponent is precisely given by the anomalous dimension 
of the Luttinger liquid that would exist for 6 = 0. This is a very important result, because 
in realistic experimental systems the interchain hopping t± can never be completely turned 
off. Our result implies that for small but finite 6 the anomalous dimension of the Luttinger 
liquid is in principle measurable, although strictly speaking the system is a Fermi liquid. 



V. COMPARISON WITH OTHER METHODS 

As mentioned in the introduction, an alternative way to approach the problem in the 
limit of weak interchain hopping is by straight or renormalization group aided perturbation 



theory in t ± |T2| , p!3| , p!7| , pO| , p2| , p3| . In contrast to our result using higher dimensional bosoniza- 
tion these methods yield a finite value for the quasi-particle weight only for values of the 
anomalous dimensions smaller than a critical value. 



Starting point in the approximation proposed by Wen [0 and elaborated on by others 



r3L2"0 , 2"2 ) |23 ] is to expand the self-energy defined in the usual way by 



G(k,u) = (88) 

in powers of the hopping. If one defines 5E(k, uS) = S(k, uj) — E(k, u;)| t and linearizes 
the dispersion in the direction of the chains near kp, one obtains for the model where only 
intrachain electron-electron interactions are taken into account 

G{KU) = gi^u^-tX^-mKuj) ' (89) 

where g(k\\,u) is the exact interacting single chain propagator. Wen |12| and later others 
|13| , |20[| then introduce the rather crude approximation to neglect <5X completely. This is 
justified by Wen by the fact that <5E formally vanishes with a higher power in t± than the 
linear term kept. For the system to be a Fermi liquid a zero frequency pole should occur. 
This can happen for arbitrarily small hopping only if g(k\\,0) diverges. Using the exact 
results for the one-dimensional Green's function, g(k\\,0) can in fact be shown to diverge 
(for the model without and with spin) for values of the anomalous dimension 7 < 1. In 
this case the quasi-particle weight in this approximation is proportional to (t±/ Ep)' y ^ 1 ~' y \ 
and therefore vanishes when 7 reaches one. This result is in agreement with the naive 
renormalization group argument in which the weak-coupling criterion for the relevance of 
t± is used for arbitrary strong coupling. This behavior of the quasi-particle weight differs 
from our bosonization result, which yields a finite value also for 7 = 1 and larger. The 
shape of the Fermi surface for 7 < 1 in Wen's approximation is determined by the equation 
g(k\\,0) = l/£j_(k_i_) and yields a modification bkpjkp compared to the flat surface of zero 
transverse hopping of order (t^/Ep) 1 ^ 1- ^ , i.e. the shape of the Fermi surface depends 
on the strength of the electron-electron interaction. Such a dependence is missing in our 
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bosonization approach. This is probably the reason for the qualitatively different results for 
the quasi-particle weight for 7 « 1. 

If one calculates the spectral function corresponding to Wen's Green's function, it shows 
various properties in the low energy regime which look rather unphysical. If one crosses the 
Fermi surface there are sharp poles even for finite distances away from it on one side, but 
there is only continuous spectral weight on the other side. We believe our bosonization result 
to produce more reliable spectra, at least for 7CI. For a direct comparison for our model 
including electron-electron interchain interaction, it would be necessary to first generalize 
Wen's approximation to this case. This can be done be replacing g(k\\,u) in Eq. (|5§D by the 
exact Green's function for zero hopping which takes into account the interchain interaction. 
For this model the approximation to neglect <5£ completely is even more serious, as #E 
vanishes with one power in t± less than in the model with intrachain electron-electron 
interaction only. This can be easily shown diagrammatically. 

Clarke and Strong [^] have argued that even in the range 1/2 < 7 < 1 the Green's func- 
tion in Wen's approximation contains unphysical non-analyticities which indicate a break- 
down of Fermi liquid behavior already at 7 = 1/2. This seems to fit well to the other 



results from their concept of "confined coherence". It will be shown elsewhere [|3[] that 
their arguments concerning the analytical properties of Wen's Green's function are not well 
justified. 

Even if one assumes it to be correct that the Fermi surface becomes flat as 7 approaches 
one (from below), the resulting state of the system cannot be a simple Luttinger liquid, 
defined as a non-trivial metallic ground state where different instabilities mutually cancel. 
Using the parquet approach it has recently been shown that two-dimensional systems 
with flat regions on opposite sides of the Fermi surface always develop some sort of instability 
towards a phase with spontaneously broken symmetry. 



VI. CONCLUSIONS 

In this work we have used the higher-dimensional bosonization approach to study the 
problem of coupled Luttinger liquids. Our main result is that strictly speaking Luttinger 
liquid behavior exists only for t± = 0. Any finite value of the interchain hopping leads 
to a finite quasi-particle residue, which we have explicitly calculated. Nevertheless, in a 
large intermediate regime of wave-vectors and frequencies the Green's function exhibits 
precisely the same scaling behavior as for t±_ = 0. Keeping in mind that experiments are 
always performed with finite resolution, the intermediate scaling regime seems to be the 
experimentally relevant one - in this sense the interchain hopping is irrelevant. 

Although our approach is non-perturbative in the sense that an infinite number of Feyn- 
man diagrams have been summed, it is only approximate. However, we have a strong non- 
perturbative argument why all terms beyond the Gaussian approximation that have been 
ignored in Eqs. (H)— (0) are negligible in the parameter regime of interest: The closed loop 
theorem p8| , p9f , which is essentially equivalent with the Ward-identity derived by Castellani, 
Di Castro and Metzner [42]|, guarantees a cancellation of the leading infrared singularities in 



the non-Gaussian terms to all orders in perturbation theory. Of course, for linearized energy 
dispersion the shape of the Fermi surface is fixed. Therefore, if the renormalization of the 
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shape of the Fermi surface by the interaction becomes relevant in the present problem, our 
conclusions are expected to be modified. However, as discussed in Sec.[V|, at least as long as 
the value of the anomalous dimension without interchain hopping is small compared with 
unity, we expect the interaction-induced modification of the shape of the Fermi surface to be 
unimportant. It should also be mentioned that only small momentum transfers have been 
taken into account in our approach, so that possible instabilities due to 2 ^-processes have 
been neglected. Thus, although we have for simplicity taken the zero-temperature limit, our 
results are implicitly restricted to temperatures T where the system is in the normal state. 
It is easy to see for T > the expressions for the long-distance behavior of the static Debye- 
Waller factor Q a (r\\v a , 0) derived above remain correct at distances small with the thermal 
de Broglie wavelength A^ = hvpjT. Beyond this length scale we find that Q a (r\\v a ,0) is 
proportional to —\r\\\/\ t h- We therefore conclude that for (Ok)^ 1 <C A^ the intermediate 
scaling regime discussed in Sec. |IVB3| exists even at finite temperature. 



ACKNOWLEDGMENTS 

This research was supported in part by the National Science Foundation under Grant No. 
PHY94-07194 during a stay of one of us (K. S.) at the ITP-workshop on Non-Fermi-Liquid 
Behavior in Solids. V. M. is grateful to the Deutsche Forschungsgemeinschaft and the NSF 
Grant No. DMR-9416906 for financial support. The work of P. K. and V. M. was partially 
supported by the ISI Foundation and EU HC&M Network ERBCHRX-CT920020. In par- 
ticular, P. K. and V. M. would like to thank Richard Hlubina for stimulating discussions 
during a workshop at Villa Gualino (Torino), which motivated us to take a closer look at 
the nesting problem. We would also like to thank Steven Strong for communications. 



23 



REFERENCES 



[1] F. D. M. Haldane, J. Phys. C 14, 2585 (1981). 

[2] C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abrahams, and A. E. Ruckenstein, 

Phys. Rev. Lett. 63, 1996 (1989). 
[3] P. W. Anderson, Science 256, 1526 (1992); P. W. Anderson, Phys. Rev. Lett. 64, 1839 

(1990); 65, 2306 (1990); 66, 3226 (1991). 
[4] S. Chakravarty and P. W. Anderson, Phys. Rev. Lett. 72, 3859 (1994); S. Chakravarty, 

A. Sudb0, P. W. Anderson, and S. Strong, Science 261, 337 (1993); L. Yin and S. 

Chakravarty, Int. J. Mod. Phys. B 10, 805 (1996). 
[5] D. Jerome, Science 252, 1509 (1991). 
[6] B. Dardel et al, Europhys. Lett. 24, 687 (1993). 
[7] M. Nakamura et al, Phys. Rev. B 49, 16191 (1994). 
[8] G.-H. Gweon at al, J. Phys.: Cond. Mat. 8, 9923 (1996). 

[9] H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990); in Jerusalem Winter School for The- 
oretical Physics, Vol. 9, Correlated Electron Systems, edited by V. J. Emery (World 
Scientific, Singapore, 1993). 
[10] L. P. Gorkov and I. E. Dzyaloshinski, Sov. Phys. JETP 40, 198 (1975). 
[11] P. A. Lee, T. M. Rice, and R. A. Klemm, Phys. Rev. B 15, 2984 (1977). 
[12] X. G. Wen, Phys. Rev. B 42, 6623 (1990). 

[13] C. Bourbonnais and L. G. Caron, Int. J. Mod. Phys. B 5, 1033 (1991). 
[14] H. J. Schulz, Int. J. Mod. Phys. B 5, 57 (1991). 

[15] M. Fabrizio, A. Parola, and E. Tosatti, Phys. Rev. B 46, 3159 (1992). 

[16] F. V. Kusmartsev, A. Luther and A. Nersesyan, JETP Lett. 55, 692 (1992); A. Ners- 

esyan, A. Luther, and F. V. Kusmartsev, Phys. Lett. A 179, 363 (1993). 
[17] C. Castellani, C. Di Castro and W. Metzner, Phys. Rev. Lett. 69, 1703 (1992). 
[18] V. M. Yakovenko, JETP Lett. 56, 5101 (1992). 
[19] A. Finkelshtein and A. I. Larkin, Phys. Rev. B 47, 10461 (1993). 
[20] D. Boies, C. Bourbonnais, and A.-M. S. Tremblay, Phys. Rev. Lett. 74, 968 (1995). 
[21] P. Kopietz, V. Meden and K. Schonhammer, Phys. Rev. Lett. 74, 2997 (1995). 
[22] D. G. Clarke and S. P. Strong, J. Phys. Cond. Mat. 8, 10089 (1996); preprint |con^ 



mat/9607141 



A. M. Tsvelik, preprint |cond-mat / 9607209| . 
[24] A. Luther, Phys. Rev. B 19, 320 (1979). 

[25] F. D. M. Haldane, Helv. Phys. Acta. 65, 152 (1992); in Proceedings of the Int. School 
of Physics "Enrico Fermi", Course 121, (North Holland, Amsterdam, 1994). 

[26] A. Houghton and J. B. Marston, Phys. Rev. B 48, 7790 (1993); A. Houghton et al, 
ibid. 50, 1351 (1994); A. H. Castro Neto and E. Fradkin, Phys. Rev. Lett. 72, 1393 
(1994); Phys. Rev. B 49, 10877 (1994); ibid. B 51, 4084 (1995). 

[27] J. Frohlich, R. Gotschmann, and P. A. Marchetti, J. Phys. A 28, 1169 (1995). 

[28] P. Kopietz and K. Schonhammer, Z. Phys. B 100, 259 (1996); P. Kopietz, J. Hermis- 
son and K. Schonhammer, Phys. Rev. B 52, 10877 (1995); P. Kopietz, Calculation of 
the single-particle Green's function of interacting fermions in arbitrary dimension via 
functional bosonization, in Proceedings of the Raymond L. Orbach Symposium, edited 
by D. Hone (World Scientific, Singapore, 1996). 



24 



[29] P. Kopietz, Bosonization of Interacting Fermions in Arbitrary Dimensions, to appear 

in Lecture Notes in Physics m: Monographs (Springer, Heidelberg, 1997). 
[30] P. Kopietz and G. E. Costilla, Phys. Rev. Lett. 76, 4777 (1996). 
[31] J. Solyom, Adv. Phys. 28, 201 (1979). 
[32] H. J. Schulz, J. Phys. C 16, 6769 (1983). 
[33] K. Penc and J. Solyom, Phys. Rev. B 47, 6273 (1993). 
[34] H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

[35] J. M. Ziman, Principles of the Theory of Solids, (Cambridge University Press, Cam- 
bridge, 1964), p.39. 
[36] A. Luther and I. Peschel, Phys. Rev. B 9, 2911 (1974). 
[37] V. Meden and K. Schonhammer, Phys. Rev. B 46, 15753 (1992). 
[38] J. Voit, Phys. Rev. B 47, 6740 (1993). 

[39] K. Schonhammer and V. Meden, Phys. Rev. B 47, 16205 (1993). 

[40] C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and 

Engineers, (McGraw-Hill, New York, 1978). 
[41] R. Hlubina, Phys. Rev. B 50, 8252 (1994). 

[42] C. Castellani, C. Di Castro and W. Metzner, Phys. Rev. Lett. 72, 316 (1994). 
[43] K. Schonhammer, to be published. 

[44] A. Zheleznyak, V. Yakovenko, and I. Dzyaloshinskii, preprint |cond-mat / 9609 1 18 . 



25 



FIGURES 





1 / 

: / 7i 

/ a i 

/ i 
' i 

i 

i \ 


K 1 


\ 


\ 






i 

T 
1 
1 

I 1 
\ 1 
\ 1 




) 





I 



FIG. 1. Fermi surface of an array of chains with small interchain hopping. Only the intersec- 
tion with the plane k z = is shown. 
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FIG. 2. Approximation of the Fermi surface in Fig.|T] by M = 8 flat patches. 
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FIG. 3. Squat boxes associated with the first two patches in Fig.||. The dashed circle indicates 
the maximal possible momentum transfer q c of the interaction. The vector k Q points to the origin 
of a local coordinate system centered at patch P£, and v Q is the corresponding local Fermi velocity. 




FIG. 4. Fermi surface of an array of chains without interchain hopping. Only the intersection 
with the plane k z = is shown. 
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FIG. 5. Fermi surface of the 4-patch model. Only the intersection with a plane of constant k z 
is shown. 
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FIG. 6. Fermi surface of the 4-patch model with curved patches. If the component of q perpen- 

q 2 

dicular to v a is denoted by q±, the patches can be described by energy dispersions £q — v a • q+ 2mj_ 
with negative effective mass m±. 
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FIG. 7. The anomalous dimension 7^ as a function of the coupling constant g = e 2 /(ttUvf) for 
two values of a±/ao (dashed and dotted lines) and / q = 4ne 2 /(aj_q) 2 (solid line). 
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FIG. 8. The Lorentzian broadened spectral function p+ for the coupled chains as a function 
ll> for different gy. The parameters are g = 0.8, a_i_/ a o = 3 and the broadening x/{ v fk) = 0.01. 
is measured relative to fi and relative to k,F- 
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